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1  SUMMARY 


This  work  evaluates  the  viability  of  an  intennediate  band  solar  cell  design,  wherein  a 
superlattice,  comprising  lattice-matched  layers  of  electronically  mismatched  GaAsl-x(Sb)Nx 
(x<0.04)  and  thin  barriers  of  AlxGal-xAs  (x<0.3),  is  inserted  within  the  intrinsic  region  of  a 
wide  bandgap  AlGaAs  p-i-n  diode.  In  the  proposed  design  the  upper  conduction  band  of 
GaAsSbN  is  maintained  in  resonance  with  the  barrier  conduction  band  of  host  material  to 
promote  an  efficient  extraction  of  electrons,  whereas  the  3D  nature  of  the  this  band  favors  a 
strong  intermediate  (superlattice  fundamental  miniband)  -  to  band  second  photon  absorption. 

The  influence  of  device  design  parameters  upon  efficiency  are  extracted  using  drift  diffusion- 
based  simulations,  that  incorporates  calculation  of  the  absorption  properties  of  the  superlattice 
(SL)  region  and  the  host  AlGaAs  p-i-n  solar  cell,  along  with  realistic  semiconductor  parameters. 
It  is  shown  that  through  a  careful  device  design  optimization  practical  1  sun  efficiencies 
exceeding  33%  (AMO),  36%  (AM  1.5)  may  be  achieved.  Analyses  of  device  performance,  under 
concentrated  sunlight,  show  the  potential  for  45%  efficiency  under  500  suns  illumination.  An 
investigation  of  the  radiation  tolerance  of  the  proposed  device  suggests  that  End  of  Life 
efficiencies  in  excess  of  29%  AMO  are  possible. 

The  study  suggest  that  this  new  class  of  highly  efficient  and  radiation  resistant,  single  junction 
solar  cells  could  deliver  practical  performance  comparable  or  exceeding  that  of  existing  multi¬ 
junctions  devices.  Considering  the  device  absorber  thicknesses  and  reduced  complexity  in  terms 
of  current  matching/lattice  matching,  the  proposed  device  also  offers  a  path  for  significant 
improvements  in  tenn  of  specific  efficiencies,  reliability  and  manufacturing  costs. 
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2  INTRODUCTION 


The  intermediate  band  solar  cell  (IBSC)  presents  a  somewhat  simpler  single-  junction  device 
alternative  to  multijunction  devices,  where  the  presence  of  an  intermediate  band  between  the 
conduction  and  valence  bands  allows  electrons  to  be  promoted  either  directly  to  the  conduction 
band  via  a  single  high-energy  photon,  or  via  the  intermediate  band  with  two  low-energy  photons 
[1].  As  only  the  population  of  two  lower  energy  transitions  (valance  band  to  IB  and  IB 
to  conduction  band)  have  to  satisfy  equality,  the  current-matching  constraint,  by 
comparison  to  multijunction  devices,  is  partially  relaxed.  The  cell  has  the  voltage  of  the 
highest  energy  bandgap,  provided  the  bands  are  thermally  isolated.  The  cell  can  also  capture 
the  solar  energy  from  low-energy  photons.  The  theoretical  efficiency  of  an  ideal  IBSC 
device  is  predicted  to  approach  45%  and  63%  under  1  sun  and  maximum  sunlight 
concentration,  respectively  [1].  Thus  far  efforts  to  produce  IBSCs  have  mainly  focused  on 
quantum  dot  (or  quantum  well)  technologies  based  on  conventional  compound 
semiconductors  (i.e.  InGaAs/GaAs  material  system)  which  have  not  yet  shown  substantial 
IBSC  (2  photons)  photo-generation  [2-4] 

Recently  strong  two  photon  absorption/IB  SC  behavior  has  been  reported  for  highly 
electronically  mismatched  bulk  alloys  [5,6]  like  in  dilute  nitrogen  alloys  of  GaAs  (GaAsN).  In 
such  alloys  the  interaction  between  the  highly  localized  impurity-like  states  of  nitrogen  with  that 
of  the  conduction  band  of  the  host  crystal  leads  to  a  splitting  of  the  conduction  band  with  an 
upper  band  denoted  here  E+  and  a  lower  branch  denoted  E-.  The  E-  conduction  band  is 
responsible  for  the  unusual  band  gap  shrinkage  observed  for  the  dilute  nitride  alloys  (—150 
meV/percent  of  nitrogen)  and  the  increase  of  electron  effective  masses  that  are  at  the  origin  of 
exceptionally  strong  near  band-edge  absorption  properties  of  these  alloys.  Unfortunately  the 
introduction  of  nitrogen  in  these  alloys  is  often  accompanied  with  severe  shortening  of  minority 
carrier  diffusion  lengths  commonly  believed  as  the  main  cause  for  open  circuit  voltage 
degradations  in  dilute  nitride  devices. 

Building  on  the  encouraging  results  shown  by  Lopez  and  co-workers  [5],  here  we  evaluate  a 
superlattice  intermediate  band  solar  cell,  wherein  a  superlattice  comprising  lattice  matched  layers 
of  electronically  mismatched  alloys  and  thin  barriers  are  inserted  within  the  intrinsic  region  of  a 
wide  bandgap  p-i-n  diode.  In  the  proposed  design  (Figure  1)  the  upper  conduction  band  of  the 
mismatched  alloys  is  maintained  to  be  in  resonance  with  the  barrier  and  bandgap  of  host 
material  to  promote  an  efficient  extraction  of  electrons,  where  the  choice  of  shallow  valence 
band  offsets  favors  the  minority  hole  extraction,  and  the  intermediate  level  is  built 
through  superlattice  minibands  formed  by  coupling  lower  band  gap  (lower  conduction 
branch  of  mismatched  alloys)  wells.  Here  carriers  can  be  promoted  either  directly  to  the 
conduction  band  or  via  the  intermediate  band,  permitting  the  absorption  of  low  energy  photons 
whilst  maintaining  a  high  cell  voltage.  Also  the  characteristic  lengths  of  the  wells  can  be 
chosen  to  be  substantially  smaller  than  typical  diffusion  lengths  of  the  electronically 
mismatched  alloy  [7],  an  approach  has  been  shown  to  overcome  the  minority  carrier  losses 
observed  in  bulk  materials  [8].  To  attain  the  necessary  combination  of  high  and  low  bandgaps 
and  low  dislocation  density,  we  use  materials  that  are  lattice-matched  to  a  GaAs  (or  Ge) 
substrate. 
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Figure  1.  Energy  diagram  (left)  and  schematic  view  (right)  of 
the  proposed  device. 


Within  the  framework  of  this  investigation  we  have  undertaken  a  modeling  of  the  superlattice 
absorption  properties  and  accounted  for  these  properties  to  derive  both  somewhat  idealist 
detailed  balance  and  practical  drift-diffusion  type  efficiency  analyses.  The  study  also  included 
gauging  the  potential  radiation  tolerance  of  the  proposed  design.  The  modeling  effort  was 
completed  by  also  evaluating  within  the  parameter  design  space  best  combination  for  maximum 
efficiency  under  1  sun  and  concentrated  sunlight  conditions. 
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3  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 


Method  assumptions  and  procedures  both  for  the  modeling  and  experimental  component  of  the 
work  are  summarized  below. 

3.1  Material  selection  and  modeling  of  the  superlattice  absorption  properties 

The  information  that  follows  in  Sections  3. 1.1 -3. 1.3  provides  useful  information  for  the  body  of 
the  technical  report. 

3.1.1  Material  selection  and  band-structure  calculations 

To  attain  the  necessary  combination  of  high  and  low  bandgaps  and  low  dislocation  density,  we 
use  materials  that  are  lattice-matched  to  a  GaAs  (or  Ge)  substrate.  In  our  design  we 
contemplate  GaAsl-x-yNxSby  or  GaAsN  quantum  well  layers  that  are  lattice  matched  to 
GaAs  and  will  be  incorporated  into  direct  bandgap  low  A1  content  (x<30%))  AlGaAs  host 
material,  to  attain  high  bandgaps  of  1.8-  1.9  eV,  and  low-energy  bandgaps  of  1-1.2  eV. 
In  fact  the  anti-crossing  interaction  between  the  N-localized  state  and  the  conduction  band  of 
host  III-V  matrix.  The  new  conduction  sub-bands  originated  by  such  a  splitting  due  to  the 
N  induced  perturbation  are  denoted  as  E_  and  E+  and  are  given  as  [9], 

Et  (k)  =  +  Em_y  (k)±  V(£„  -  Em_f  (k)f  +  4xV„2  }  (1) 

AlGaAs  GaAsN 


Percentage  of  Aluminum  Percentage  of  Nitrogen 


Figure  2.  Calculated  band  structure  for  AlGaAs  &  GaAsN 
material  system  &  example  of  composition  selection  A1  and  N. 

where  En  is  the  energy  of  the  localized  nitrogen  state,  EIII-V(k)  is  the  dispersion  of  the  host 
crystal  conduction  band,  Vn  is  the  strength  of  the  anti-crossing  interaction  between  the  nitrogen 
localized  states  and  the  conduction  band  states  of  host  matrix,  and  x  is  the  nitrogen 
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concentration.  All  the  energies  are  relative  to  the  valence  band  maximum  of  the  unperturbed 
host  crystal.  The  band  anticrossing  in  GalnNAs  Alloys  (BAC)  model  with  VN=2.7eV  provides 
an  excellent  fit  for  the  nitrogen  concentration  x~0.04.  The  E-  and  E+  conduction  band 
branches  of  the  dilute  nitride  alloy  have  an  energy  separation  that  is  generally  set  by  the 
nitrogen  composition  of  the  dilute  nitride  material.  The  calculated  band  structure  for 
AlGaAs  and  GaAsN  material  system  is  shown  in  Figure  2.  The  shaded  minibands  in  Figure 
2  correspond  to  fundamental  electron  (green)  and  heavy  hole  (gray)  bands  in  3  mn  GaAsN 
(N'2.5%). 

The  use  of  thin  barriers  and  wells  can  be  used  to  adjust  the  miniband  energies  hence  providing 
an  additional  degree  of  freedom  in  the  energy  of  the  intermediate  band  energy  independently 
from  the  E+  energy  level.  Here  the  shallow  valence  band  offsets  favor  a  rapid  thermionic 
escape  of  holes  at  a  time  scale  of  a  few  picoseconds.  Note  considering  the  deep  electron  wells 
the  thenno-ionic  escape  of  electrons  from  the  intermediate  level  exceeds  the  nanosecond  scale 
and  hence  the  only  possible  escape  route  is  the  absorption  of  the  second  photon. 

3.1.2  Calculation  of  the  absorption  properties  of  the  superlattices. 

The  calculation  of  the  SL  region  absorption  properties  follow  a  previously  established 
methodology  [10,  11].  The  methodology  includes  the  extraction  of  the  superlattice  miniband 
energies  using  the  transfer  matrix  method  and  a  formalism  detailed  in  Appendix  A  is  then  applied 
to  extract  the  absorption  coefficient.  The  formalism  associated  with  calculation  of  miniband 
energies  and  the  absorption  coefficient  is  detailed  in  Appendix  A. 


Figure  3.  Absorption  coefficient  extracted  from  spectroscopic 
ellipsometry  measurements  of  a  pseudomorphically  strained 
bulk  GaAsN  grown  on  GaAs  (100). 

The  simulation  effort  included  modification  of  the  existing  Matlab-based  custom  source  code  for 
superlattice  to  include  material  parameters  for  dilute  nitride  and  AlGaAs  alloys.  The  calculation 
account  for  strain  Hamiltonian  and  band-anti-crossing  effects.  Absorption  of  the  bulk  dilute 
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nitride  was  extracted  from  spectroscopic  ellipsometry  measurements  of  previously  grown  bulk 
pseudomorphically  strained  layers  of  1.2  eV  GaAsN  shown  in  Figure  3. 

An  example  of  the  calculation  of  the  absorption  properties  of  the  intennediate  band  superlattice 
structure  incorporating  30  period  of  a  5nm  thick  GaAsN  wells  and  5nm  thick  AlGaAs  barriers  is 
shown  in  Figure  4. 


Energy  (eV) 


Figure  4.  Calculated  absorption  coefficient  of  a  30  period 
GaAs0.98N0.02  (3nm)/  Alo.20Gao.80As  (3nm)  Superlattice  following 
the  methodology  developed  in  ref  10. 

3.2  Simulation  of  Device  Performance. 

This  portion  of  the  modeling  effort  relied  on  simulations  of  the  device  I-V  characteristics  and 
efficiencies  using  a  Detailed-Balance  (thermodynamic  somewhat  idealist  model)  and  Drift 
Diffusion  based  approach.  Existing  Matlab-based  custom  codes  were  modified  to  incorporate 
the  chosen  parameters  and  assumption  detailed  hereafter  in  the  code 

3.2.1  Detailed  Balance  Simulations 

The  detailed  balance  modeling  approach  was  similar  to  the  one  described  in  reference  l,  except 
that  to  add  a  higher  degree  of  realism  the  modeling  accounted  for  the  absorption  properties  and 
finite  thicknesses  of  SL  intermediate  band  (IB)  region.  To  simulate  the  efficiencies  we  have 
extrapolated  absorption  properties  calculated  in  Figure  4.  These  calculations  were  performed 
both  under  space  sunlight  illumination  (AMO)  and  terrestrial  (AM1.5G)  conditions. 
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3.2.2  Drift-Diffusion  simulations  and  modeling  of  radiation  effects 

In  order  to  gain  a  more  realistic  insight  on  the  potential  of  the  proposed  design  we  have 
undertaken  a  drift-diffusion  modeling  as  a  function  emitter  and  base  thicknesses  and  doping 
and  by  incorporating  realistic  minority  carrier  lifetimes  and  Shockley  Reed  Hall  (SRH) 
recombinations  using  a  previously  established  methodology  [12].  For  the  superlattice  region 
we  have  assumed  a  superlattice  comprising  N  period  (N  varied  from  10-100)  of  5nm  AlGaAs 
barriers  and  5nm  GaAsN(Sb)  wells.  The  simulations  were  carried  out  with  AlGaAs  n-type  base 
doping  range  of  10  to  10  cm"  while  the  p-type  emitter  doping  range  is  between  10  and 
5*1018cm"3.  SRH  lifetime  for  depleted  region  is  between  10'8  sec  to  5*10"8  sec  while  the 
surface  recombination  velocity  is  in  the  order  of  103cm/sec.  For  uniformity  in  results  it  is 
assumed  3%  anti-reflective  coating  and  5%  grid  shadowing  losses  have  been  included. 
Minority  carrier  lifetimes  as  a  function  of  doping  concentration  for  GaAs  instead  of 
AlGaAs  are  taken  including  Shockley-Read-Hall  effect,  band  to  band  recombination  and  Auger 
recombination.  Minority  carrier  diffusion  lengths  and  diffusion  coefficient  as  a  function  of 
doping  concentrations  and  absorption  coefficient  of  material  are  taken  from  the  experimental 
values  in  the  literature. 

Here  we  also  attempted  to  simulate  the  effect  of  space  radiation  induced  degradations  by 
assuming  a  radiation  induced  degradation  of  minority  carrier  diffusion  length  similar  to  that  of 
GaAs.  The  reduction  in  minority  carrier  diffusion  length  can  be  modeled  according  to  relation, 

L"2=Lo"2+Kl  (p  (2) 

Where  L  is  the  net  diffusion  length,  Lo  is  the  diffusion  length  of  defect  free  III-V  material,  cp  is 
the  electron  radiation  fluence  per  cm"  and  KL  is  radiation  constant,  which  depends  on  the 
material  used  in  the  sub-cell.  The  Kl  is  taken  equal  to  3x10-8  for  fluencies  of  1  MeV  electrons.  It 
should  be  noted  that  analyses  presented  here  do  not  account  for  carrier  removal  effects  which  in 
turn  can  affect  the  E-field  across  the  superlattice  region. 
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4  RESULTS  AND  DISCUSSION 


4.1  Detailed  balance  simulation  of  the  device  accounting  for  SL  absorption 

To  assess  the  photovoltaic  potential  of  the  proposed  design  we  have  undertaken  a  detailed 
balance  modeling  by  assuming  a  device  incorporating  a  0.5  microns  superlattice  region  with 
the  intermediate  band  region  absorption  properties  comparable  to  the  one  calculated  in  Fig  3. 
The  preliminary  calculation  presented  here  assumes  that  all  photons  with  energy  larger 
than  the  bandgap  of  AlGaAs  (or  E+)  are  absorbed  in  the  emitter  of  the  device  (thick 
emitter  approximation)  whereas  photons  with  energy  comprised  between  the  fundamental 
miniband  state  and  resonant  E+  level  are  absorbed  with  absorption  properties  similar  to  that 
of  the  dilute  nitride  superlattice.  The  absorption  strength  for  inter-subband  photon 
absorptions  between  the  fundamental  and  excited  states  of  the  superlattice  has  not  been 
considered  to  account  for  the  fact  that  under  normal  light  incidence  the  transition  from 
fundamental  to  first  excited  miniband  state  is  forbidden  and  the  oscillator  strength  for  such 
transition  in  a  pure  2D  state  is  zero.  Since  here  the  inter-subband  absorption  between  the 
fundamental  or  excited  states  (intermediate  2D  band)  to  the  upper  conduction  band  and 
involves  3D  like  E+  states  of  the  dilute  nitride  material/AlGaAs,  thus  the  oscillator 
strength  =0  selection  rule  does  not  limit  the  absorption  as  is  the  case  for  conventional 
superlattices.  For  the  calculations  presented  here  we  have  assumed,  that  within  the  limit  of 
available  photons  for  each  photon  absorbed  by  the  IB  a  second  photon  is  absorbed  between  IB 
and  E+. 


Figure  5.  Detail  balance  theoretical  efficiency  of  an  IBSC 
device  with  inclusion  of  absorption  properties  of  a  0.5  micron 

SL  region. 

The  AMO  efficiency  contour  plots  as  a  function  of  the  superlattice  fundamental  minband 
energy  and  the  host  bandgap  for  a  superlattice  IBSC  device  proposed  here  are  shown  in  figure  5 
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(calculation  accounts  for  10%  reflection  and  grid  shadowing  losses).  The  dotted  rectangle  in 
Figure  5  is  for  eye  guidance  and  indicates  the  practical  design  space  available  with  the  proposed 
AlGaAs/GaAsSbN  SL  material  system  described  in  section  3-1.1.  These  somewhat  idealist 
calculations  stress  the  potential  of  the  proposed  IBSC  for  achieving  1  sun- AMO  efficiencies  in 
excess  of  39%. 


Figure  6.  Detailed  Balance  efficiency  contour  plots  for  1000 
sun  illumination  as  a  function  of  Superlattice  energy  (Ei)  and 
bandgap  of  host  cell  material  (Eg) 


Next  we  have  extended  the  efficiency  calculations  under  concentrated  sunlight.  In  Figure  6  the 
efficiency  contour  plots  for  the  proposed  IBSC  device  are  computed  assuming  a  0.5  micron 
thick  superlattice  region.  This  detailed  balance  simulation  suggest  sets  an  upper  theoretical 
limit  exceeding  55%  for  a  device  operating  under  1000  sun  AM  1.5  illumination. 

4.2  Drift  Diffusion  Efficiency  Simulations 

In  order  to  gain  a  more  realistic  assessment  of  the  potential  of  the  proposed  design  we  have 
undertaken  the  development  of  a  drift-diffusion  model  as  a  function  of  emitter,  base  thickness, 
doping,  and  incorporating  realistic  minority  carrier  lifetimes  and  SRH  . 

4.2.1  Optimizing  superlattice  and  emitter  thicknesses 

For  the  superlattice  region  we  have  assumed  a  superlattice  comprising  N  period  (N  varied  from 
10-200)  of  5nm  AlGaAs  barriers  and  5nm  GaAsNSb  wells.  Figure  7  shows  the  evolution  of  the 
AMO  efficiencies  as  a  function  of  the  i-region  thickness  and  emitter  thickness  with  a  fixed  base 
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thickness  of  0.5  microns.  For  a  1.8  eV  AlGaAs  an  80-90  period  SL  with  an  emitter  thickness  of 
0.5micron  provides  the  highest  AMO  efficiency  of  3 1 .4%. 


Figure  7.  Drift  Diffusion  simulation  of  the  proposed  IB  SC  for  a 
1.8eV  AlGaAs  with  1.2  eV  GaAsN/AlGaAs  superlattice  region. 

It  should  be  noted  that  due  to  strong  below  AlGaAs  bandgap  absorption  in  the  superlattice 
region,  in  the  device  simulation  for  superlattice  exceeding  a  thickness  of  0.2  microns  the  device 
IV  characteristics  were  found  to  be  insensitive  to  the  device  base  thickness.  The  IV 
characteristics  for  the  best  device  in  these  conditions  are  shown  in  an  inset  in  Figure  7. 

4.2.2  Optimizing  the  bandgaps  of  the  superlattice  and  the  AlGaAs  host  cell 

Next  we  have  evaluated  the  effect  of  the  miniband  energy  and  the  bandgap  of  host  cell  upon 
the  evolution  of  device  efficiency.  A  set  of  contour  plots  were  generated  as  a  function  of  the 
bandgap  of  the  superlattice  and  the  AlGaAs  host  bandgap.  For  the  superlattice  region 
individual  layer  thicknesses  were  adjusted  to  yield  the  targeted  minband  energies  of  1.1,  1.2  and 
1.3  eV.  Figure  8  represents  the  results  of  these  calculations.  Each  data  point  in  the  Figure 
corresponds  to  the  maximum  efficiency  extracted  from  full  contour  plot  calculation  similar  to 
those  shown  in  Fig  7. 

Interestingly  it  is  found  that  the  AMO  efficiency  can  be  further  boosted  by  using  a  1.9 
eV  AlGaAs  material,  with  a  targeted  efficiency  of  -34%  AMO  (Figure  8a).  It  should  be 
noted  that  from  a  practical  standpoint  it  may  be  preferable  to  target  a  slightly  lower  host 
bandgap  i.e.  1.85  eV  which  could  yield  an  efficiency  of  -33%  (Figure  8b). 


Approved  for  Public  Release;  distribution  is  unlimited. 

10 


Figure  8.  (a)  Evolution  of  AMO  Efficiency  as  a  function  of 
AlGaAs  bandgap.  (b)  a  contour  plot  with  AlGaAs  bandgap  of 
1.85eV  and  miniband  energy  at  1.2  eV. 

4.2.3  Effect  of  Sunlight  Concentration  (terrestrial  applications) 

We  have  also  evaluated  the  potential  of  the  technology  for  terrestrial  and  concentrator 
applications.  Figure  9  shows  the  best  efficiency  recorded  for  AM1.5G  illumination.  A  1  sun 
efficiency  of  ~36%  seems  achievable  by  using  a  1.2  eV  superlattice  and  1.85  eV  AlGaAs.  In 
Figure  9,  each  period  of  the  superlattice  is  lOnm  (i.e.  0.1microns=  10  periods). 
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Figure  9.  Drift-Diffusion  Simulation  of  the  AMI. 5  efficiency 
as  function  of  concentration  and  SL  region  thickness. 

It  is  worth  noting  that  these  max  efficiencies  are  reached  for  i-region  thickness  of  ~  2  microns 
(100  periods)  assuming  a  fully  depleted  i-region  which  assumes  a  background  carrier 
concentration/unintentional  doping  in  the  I  region  lower  or  about  ~10  cm'  ). 

4.2.4  Evaluation  of  Device  End-of  Life  Efficiencies 

We  next  expanded  our  simulations  to  include  the  effect  of  radiation  degradation  to  assess  the  end 
of  life  (EOL)  efficiencies  of  these  devices  in  space.  Figure  10  represents  the  evolution  of 
maximum  achievable  efficiency  for  variation  radiation  fluencies.  Each  calculation  point  is  the 
outcome  of  the  optimization  of  all  design  parameters.  The  dotted  line  in  Figure  10  indicates  the 
data  point  extracted  from  the  contour  plot  and  corresponding  IV  shown  in  Figure  1 1 . 
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Figure  10.  Evolution  of  max  AMO  end  oflife  efficiency  as  a 
function  of  IMeV  equivalent  radiation  fluence  doses. 

As  an  example  for  a  radiation  dose  of  2x1015  cm-2,  that  corresponds  to  a  degradation  similar  to 
10  years  operation  in  a  geostationary  orbit,  a  contour  plot  as  a  function  of  emitter  and 
superlattice  region  for  1.2  eV  superlattice/1.8eV  AlGaAs,  is  shown  in  Figure  11.  For  the  latest 
radiation  dose  the  study  indicates  the  potential  to  exceed  EOL  efficiencies  of  29%  AMO. 
It  should  be  noted  that  the  calculations  do  not  account  for  carrier  removal  effects  which  could 
lead  to  a  radiation  induced  increase  of  background  carrier  concentrations  in  the  i-region  which 
in  turn  may  affect  the  carrier  extraction  processes  from  the  i-region  (E-field) 
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Figure  11.  (a)  Contour  plot  of  device  EOL  efficiency  as  a  function  of  emitter 
and  i-region  thickness  for  a  IMeV  electron  fluence  dose  of  2xl015cm-2  (b)  EOL 
I-V  characteristic  of  the  device. 

4.3  Discussion  on  Effect  of  relative  absorption  strength  for  inter-sub-band  transitions 

The  basic  operation  of  an  IBSC  solar  cell  requires  a  strong  absorption  from  the  valence  band  to 
the  intermediate  band  IB,  a  condition  satisfied  by  the  superlattice  structure  (Fig.  4),  followed  by 
an  equally  strong  aborption/photo-promotion  of  carriers  (here  electrons,  as  hole  band  offset  is 
negligible)  from  the  superlattice  fundamental  miniband  to  the  host  cell  continuum.  In  practical 
realization  involving  quantum  dots  or  conventional  Kane-like  2D  superlattices  and  quantum 
wells  this  second  aborption  has  been  extremely  weak  [2-4].  The  design  contemplated  in  this 
work  builds  on  a  transition  between  the  1 st  miniband  (E-  band  like  behavior)  and  the  continuum 
that  is  resonant  with  3D  like  E+  states  populated,  and  hence  as  shown  for  bulk  like  dilute 
nitrides,  the  2nd  photon  induced  electron  promotions  are  stronger  [5].  Yet  it  is  important  to 
gauge  the  actual  strength  of  such  inter-subband  carrier  promotion,  as  it  requires  both  for  the 
band  to  be  populated  and  for  the  associated  absorption  coefficient  to  be  strong.  In  the 
calculations  presented  in  section  4.2  we  have  assumed  that  the  absorption  strength  was  similar 
to  that  of  the  miniband  which  may  be  a  reasonable  assumption  as  E-  and  E+  states  have 
somewhat  similar  density  of  states,  yet  the  overlap  of  wave  functions  (oscillator  Strength)  for  the 
transition  from  E-  or  SL  miniband  to  the  upper  conduction  band  E+  could  be  much  weaker 
than  that  of  the  inter-band  one  (valence  band  (VB)-to  -superlattice  miniband).  In  the  absence 
of  actual  measures  for  intersubband  absorption  coefficients  in  Figure  12  we  represent  the 
effect  of  the  absorption  coefficient  mismatch  between  the  two  transitions  on  the  device 
efficiency.  It  should  be  noted  that  while  a  small  mismatch  (i.e.  factor  of  2)  has  almost  no 
effect  on  the  maximum  projected  efficiency,  for  higher  absorption  mismatches  the  efficiency 
steeply  decreases,  to  eventually  drop  below  that  of  a  single  junction  1.8  eV  solar  cells  (-20% 
AMO). 
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Figure  12.  Evolution  of  device  efficiency  as  a  function  of  the  ratio  of  absorption 
coefficients  for  transitions  from  the  SL-miniband  to  conduction  band  of  AlGaAs 
and  from  valence  band  to  SL-minband. 

The  calculation  however  assumes  a  single-pass  normal  light-path  through  the  device  and  thus 
simply  incorporating  a  reflective  mirror  at  the  back  of  the  device  should  improve  the  2  photon 
absorption/up-conversion  processes.  Ultimately,  incorporation  of  a  Lambertian  reflector  which 
could  enhance  the  effective  light  path  by  4n“  (n  index  of  refraction  ~3.4),  i.e.  40  fold  could 
be  used  to  mitigate/or  completely  overcome  the  risk  associated  with  weak  inter-sub-band 
absorption  for  these  devices. 
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5  CONCLUSIONS 


In  this  work  we  investigated  the  potential  for  space  and  terrestrial  applications  of  a  novel  single 
junction  solar  cell  where  a  superlattice  comprising  well  of  1-1.3  eV  dilute  nitride  alloys  and  (1.7- 
1 .9eV)  AlGaAs  barriers  is  inserted  in  the  intrinsic  (i)  region  of  an  AlGaAs  p-i-n  diode.  The  work 
included  the  simulation  of  band-structure  and  absorption  properties  of  the  superlattice.  The  latest 
were  incorporated  in  a  Detailed  Balance  and  Drift-Diffusion  simulator  to  extract  somewhat 
ideal  and  practical  efficiency  limits  of  the  device.  The  study  included  exploring  a  wide  range  of 
design  parameters 

The  analyses  indicated  the  potential  for  this  single  junction  device  to  deliver  efficiencies  in  excess 
of  33%  under  space  sunlight  illumination  (AMO)  and  46%  under  1000  sun  concentration.  The 
effects  of  radiation-induced  degradations  were  also  incorporated.  The  potential  for  the  device  to 
deliver  end-of-life  efficiencies  in  excess  of  29% AMO  for  a  radiation  fluence  of  10  cm'~  IMeV 
electron  was  predicted. 

It  should  be  noted  that  the  technology  provides  a  path  to  achieve  efficiencies  comparable  to 
4  junction  device  in  several  fold  thinner  single  junction.  The  design  also  relieves  to  some  extent 
the  current-matching  constrains  of  multijunction  devices.  Consequently  by  comparison  to 
current  state-of-the-art,  the  technology  has  the  potential  of  enhancing  both  specific 
efficiencies  and  reliability  of  high  efficiency  space  solar  cells. 
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6  RECOMMENDATIONS 


The  device  design  studied  here  appears  as  promising  however  a  key  parameter  toward  the 
success  of  the  approach  is  to  the  absorption  strength  of  the  intersubband  transition  from  the 
superlattice  miniband  to  the  continuum.  In  our  preliminary  modeling  we  have  shown  the 
criticality  of  the  interband  absorption  toward  exceeding  the  efficiency  limit  of  single 
junction  devices.  One  approach  to  further  strengthen  the  inter-subband  aborption  could 
consist  of  enhancing  the  optical  path  of  light  by  incorporating  a  back  reflector  (i.e.  gold)  that 
recycles  the  transparency  losses  into  the  superlattice  region.  It  should  be  noted  that 
incorporation  of  a  Lambertian  reflector  could  potentially  result  in  a  4  n“  enhancement  of  the 
effective  absorption  path.  For  III-V  compounds  contemplated  here  (index  of  refraction  n~3.4) 
this  could  present  a  40  fold  enhancement  of  the  absorption  strength.  Development  of  such  a 
reflector  will  also  allow  the  device  superlattice  region  to  be  made  much  thinner  and  hence 
in  principle  will  improve  the  device  radiation  to lerance/end  of  life  efficiencies. 
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APPENDIX:  METHODOLOGY  FOR  CALCULATON  OF  MINBAND 
ENERGIES  AND  ABSORPTION  COEFFICIENT  OF  A  SUPERLATTICE 


A.l  Miniband  Energy  Calculation 


For  a  fixed  energy  and  fixed  in-plane  wave  vectors  kx  and  ky,  we  can  construct  the  total  envelope 
function  of  a  superlattice  as  a  linear  combination  of  the  eigenvectors  of  the  corresponding 
Hamiltonian.  For  each  region  j,  the  full  envelope  function  has  the  general  shape 
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where  the  sum  is  over  all  the  bands  q,  conduction  band  and  valence  band,  a/  and  bj  are  complex 
constants,  Fq'  are  eigenvector  matrices  and  /,■  defines  the  jlh  interface  from  the  origin. 
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Boundary  conditions  at  each  interface  are 


\p./  _  vp  j+ 1 

Z=l  :  Z=l 


and 


m , 


dxVj 

+ 

. 

1 

dz 

mi«  & 

(A2) 


Z=l  : 


Here,  in,  is  the  effective  mass  in  jth  layer.  For  kx=ky=  0,  the  mj=± 3/2  hole  band  gets  decoupled  from 
the  rest  of  the  bands,  and  transfer  matrix  can  be  found  applying  above  two  boundary  conditions 
to  relate  constants  a  and  b  of  the  jth  and  (j+1) th  layers,  and  is  given  as 


where 


\PJJ 


-  Mj(j+ 1) 


fa^ 

\bJ+1  J 


(A3) 


M 


J(M) 


1  +  - 


m  ,k 


j,vj+ 1 
mJ+\kj 


v 


1-- 


m  k 


j'vj+ 1 

mj+\kj 


x  mjkj+ 1  l  +  mjkJ+ 1 


mj+\kj 


mi+\kj  j 


(A4) 


Here,  Lj  is  the  thickness  of  the  jth  layer.  For  large  number  of  layers,  the  total  transfer  matrix  can 
be  written  as  the  product  of  transfer  matrices  across  each  layer 


Approved  for  Public  Release;  distribution  is  unlimited 

19 


/  l\  (  N\ 

Cl  Cl 

7  i  1  total  7  n 

Vb  )  \P  ) 
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(A5) 


We  now  require  a  =  bN  =0  for  decaying  solutions  in  the  first  and  last  barriers,  width  of  which 
is  assumed  to  be  large  enough  for  this  condition  to  hold.  To  satisfy  this  condition,  the  element 
mn  goes  to  zero  and  hence  the  solution  of  equation  mn(E)= 0  provides  the  energy  states  of  the 
corresponding  carriers  in  mj=± 3/2  hole  states.  Same  process  can  be  repeated  for  mj=±\!2  states 
also,  but  due  to  the  interaction  with  spin  orbit  split-off  states  we  have  to  deal  with  4x4  transfer 
matrix,  with  added  numerical  complexity. 

A.2  Calculation  of  the  absorption  of  the  superlattice 

Optical  absorption  in  a  semiconductor  material  is  one  of  the  very  important  physical  phenomena 
to  be  considered  for  the  application  in  many  optoelectronic  devices.  The  parameter  that 
characterizes  the  optical  absorption  is  called  the  absorption  coefficient  and  is  defined  as  the  ratio 
of  the  number  of  photons  absorbed  per  unit  volume  per  second  to  the  number  of  photons  injected 
per  unit  area  per  second.  Hence  the  absorption  coefficient  gives  the  fraction  of  photons  absorbed 
per  unit  distance  and  is  denoted  by  a  (cm'1). 

The  absorption  coefficient,  a  in  the  crystal  for  a  bulk  semiconductor  is  given  as 

~E.)  (A6) 


where,  A0  = 


2  2 

ne  i  ■  pc 


nrc£0m0a> 


i  is  the  unit  vector  along  the  direction  of  the  propagation  of  the 


incoming  photon,  pcv  is  the  transition  matrix  element  between  the  conduction  and  valence  band,  e 
is  an  electronic  charge,  nr  is  refractive  index,  c  is  the  velocity  of  the  electromagnetic  wave  in 
vacuum,  so  is  the  pennittivity  of  the  vacuum,  mo  is  the  electron  mass,  co  is  the  angular  frequency 

l  jU  c 

of  the  incoming  electromagnetic  wave,  h  is  the  Planck’s  constant/27T,  Ry  =- - — — , 

(4  ne,)  2  h 


a{)  =  — - - ,  jiir  is  the  reduced  effective  mass  of  holes  and  electrons  and  is  given  as, 

e  Mr 

—  =  —I—  +  — I— ,  £s=Kso,  K  is  a  dielectric  constant,  0„(O)  is  a  normalized  wave  function  for  both 
Mr  me  mh 

bound  and  continuum  states  and  is  given  as 
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(A7) 


bound 


3  3 

na^n 


and 


4(o)|, 


tt/Je/R 

e1^ 

4 7:Ryal  sinh (tt  JE/Ry  ) 


(A8) 


respectively.  With  inclusion  of  this,  the  total  absorption  (also  called  Elliott  formula)  due  to  both 
bound  and  continuum  states  in  bulk  semiconductor  can  be  given  as 


A 


o 


2  n2Rya\ 


(A9) 


where 


fico-  Eg 

Ry 


and  Ais  Lorenzian  width  corresponding  to  the  phonon  broadening. 


For  the  superlattice  structure  discussed  in  the  section  IV,  we  will  use  an  anisotropic  medium 
approach,  in  which  the  absorption  coefficient  of  the  bulk  material  is  modified  to  take  care  of  the 
anisotropy  introduced  by  the  presence  of  the  superlattice  periodicity.  As  discussed  in  the  section 
IV,  due  to  the  penetration  of  the  carrier  wavefunction  into  the  adjacent  wells  separated  by  thin 
barriers  in  a  superlattice,  the  coupling  between  the  two  neighboring  wells  is  very  large  to  allow 
for  the  hybridization  of  the  original  discrete  energy  levels  of  isolated  wells  into  minibands.  The 
electron  and  hole  tunneling  through  the  minibands  leads  to  the  motion  in  the  growth  direction, 
which  is  characterized  by  the  effective  masses  modified  by  the  superlattice  potential.  Superlattice 
is  then  characterized  by  a  three  dimensional  effective  medium  with  different  in-plane  and 
growth-direction  effective  masses. 

The  superlattice  absorption  coefficient  then  can  be  assumed  as  a  deviation  from  a  bulk 
absorption  coefficient  by  considering  the  anisotropy  introduced  in  the  reduced  effective  masses 
of  the  carriers  in  the  plane  of  the  layers  (in-plane  reduced  effective  mass,  ]U\\)  and  perpendicular 
to  the  plane  of  layers  (growth-direction  or  transverse  reduced  effective  mass,  respectively. 
The  ratio  of  these  two  masses  can  be  defined  as  an  anisotropy  parameter, 


7  =  /V/A  (A  10) 

This  parameter  modifies  the  absorption  coefficient  as  suggested  in  Ref.  53,  and  is  given  as 


a 


(Tico,y)  = 


A, n 


2x2r[,  (rK(/)3 


f  2  ] 

2 

l1  +  rJ 

*1 


T/n 


n=\  (*  M  +  l/«  )+r  sinh  nj -y/s'(/) 


(All) 
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where 


S'M= 


h(o-  Eg 

Kir)  ' 


Ansh 
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W =^,K(rh 


em 


1  m^eA  2  H\ 

nh  = 


(4  ns^)  2  ti“ 


\  +  y 


The  in-plane  and  transverse  reduced  effective  mass  can  be  calculated  as 


— +  —  — 


— +  — 


M\\  me  mh,  Ml  m 


(A12) 


The  in-plane  effective  mass,  m\\  is  assumed  to  be  average  of  effective  masses  of  two  different 
layers,  w  and  b  of  thickness  Lw  and  L/,  respectively,  in  superlattice 


777  W  T  -1-  777  ^  f 

m  =  — — —  with  i=e,  Ih  or  hh. 
11  I  + 1 


(A  13) 


As  a  first  approximation  transverse  effective  mass,  m±  can  be  obtain  from 


1  1  d2E 


mL  h2  dq2 


q=0 


(A  14) 


where,  energy  E  is  extracted  from  a  dispersion  relation  for  a  superlattice 

f  0 

r  d 


if  0 

cos(qd)  =  cos(kwLw)cosh(kbLb )  —  —  sinh(^Zjsin(A;wZj 


(A  15) 


here,  kw  = 


2  mwE  k  =  \2mh{V-E) 


.2  ’  b 


for  E<V  and £  = 


kJnb 
Km  H 


Using  equations  (15)  and  (16),  we  may  derive  the  transverse  effective  mass  as 


m  =  [-  cos(lwkw )sinh(Lbkb )LhB  -  sm{Lwkw)cos\i{Lbkb  )LwA ] 


tf_ 

d2 

Pr_ 

l2 


1  1  mh 


2  m 


A_+Bk± 

ku  k2; 


sin 


\\{Lbkb)sm{LwkM) 


w  Vlb  ^b  J 


(A  16) 


-  -  jjr  cosh  (Lbkb  )sm(Lwkw)LhB  +  smh(Lbkh  )cos(Lwkw)LwA } 
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where,  A  = 


B  = 


m. 


2  h2(V-E) 


and  d=Lw+Lb. 


Since  the  dielectric  constant  of  GaAs(N)  and  AlGaAs  are  not  very  different  from  each  other,  the 
ratio  of  the  in-plane  and  transverse  dielectric  constants  for  a  superlattice  structure  is  ~1.  Hence, 
the  dielectric  constant  of  the  superlattice  is  taken  as  the  weighted  average  for  two  layers  of 
superlattice  as 


£2  _  £wLw  +  £b^b 

Lw/£w+Lb/£b 


(A  17) 


For  y=l,  equation  (All)  becomes  similar  to  the  equation  (A9)  and  the  condition  corresponds  to 
the  one  of  no-anisotropy  or  the  bulk  material  absorption  coefficient. 
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LIST  OF  ACRONYMS,  ABBREVIATIONS,  AND  SYMBOLS 


Acronym/  Abbreviation  Description 


AMO 

Air  Mass  Zero 

BAC 

Band  Anti-Crossing 

EOL 

End  of  life 

IBSC 

Intermediate  Band  Solar  cell 

SL 

Superlattice 

SRH 

Shockley  Reed  Hall 

Symbol 

Description 

a 

Absorption  coefficient 

E- 

Lower  branch  of  conduction  band  of  dilute  nitride 

E+ 

Upper  branch  of  the  conduction  band  of  dilute  nitride 

Eg 

Bandgap 

<l> 

Radiation  flux 

h 

efficiency 

L 

Minority  Carrier  Diffusion  length 

KL 

Radiation  coefficient 
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